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Abstract. This work examines various statistical distributions in connection with random Vander- 
monde matrices and their extension to d-dimensional phase distributions. Upper and lower bound 
asymptotics for the maximum singular value are found to be 0(logN d ) and 0(log N d / log log N d ) 
respectively where N is the dimension of the matrix, generalizing the results in [5T]. We further 
study the behavior of the minimum singular value of a random Vandermonde matrix. In particular, 
we prove that the minimum singular value is at most N 2 exp(— Cy/N)) where N is the dimension of 
the matrix and C is a constant. Furthermore, the value of the constant C is determined explicitly. 
The main result is obtained in two different ways. One approach uses techniques from stochastic 
processes and in particular, a construction related with the Brownian bridge. The other one is a 
more direct analytical approach involving combinatorics and complex analysis. As a consequence, 
we obtain a lower bound for the maximum absolute value of a random complex polynomial on the 
unit circle. We believe that this has independent mathematical interest. Lastly, for each sequence 
of positive integers {k p }^Li we present a generalized version of the previously discussed random 
Vandermonde matrices. The classical random Vandermonde matrix corresponds to the sequence 
k p — p—1. We find a combinatorial formula for their moments and we show that the limit eigenvalue 
distribution converges to a probability measure supported on [0, oo). Finally, we show that for the 
sequence k p = 2 P the limit eigenvalue distribution is the famous Marchenko-Pastur distribution. 



1. Introduction 



Large dimensional random matrices are of much interest in statistics, where they play a roll in 
multivariate analysis. In his seminal paper, Wigner [23] proved that the spectral measure of a 
wide class of symmetric random matrices of dimension N converges, as ./V — > oo to the semicircle 
law. Much work has since been done on related random matrix ensembles, either composed of 
(nearly) independent entries, or drawn according to weighted Haar measures on classical (e.g., 
orthogonal, unitary, simplectic) groups. The limiting behavior of the spectrum of such matrices 
and their compositions is of considerable interest for mathematical physics and information theory. 
In addition, such random matrices play an important role in operator algebra studies initiated 
by Voiculescu, known now as free (non-commutative) probability theory (see [22J and [23J and the 
many references therein) . The study of large random matrices is also related to interesting questions 
of combinatorics, geometry, algebra and number theory. More recently, the study of large random 
matrices ensembles with additional structure have been considered. For instance, the properties of 
the spectral measures of random Hankel, Markov and Toeplitz matrices with independent entries 
have been studied in [6]. 

In this paper we study several aspects of random Vandermonde matrices with unit magnitude 
complex entries and their generalizations. An N x L matrix V with unit complex entries is a 
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Vandermonde matrix if there exist values 9\, . . . ,9l €. [0, 1] such that 

(\ ... i \ 



N 



= 27ri0i g27ri6i 



y e 27ri(JV-l)6»i e 27vi(N-l)e L J 

(see [17] or |21j for more details). A random Vandermonde matrix is produced if the entries of the 
phase vector 9 := . . . , 9l) G [0, 1] L are random variables. For the purposes of this paper it will 
be assumed that the phase vector has i.i.d. components, with distribution drawn according to a 
measure v which has a density fix) on [0, 1]. 

Vandermonde matrices were defined in 1131 and were also called d-fold Vandermonde matrices. The 



case d = 1 are the matrices in (1.1). For d > 2, these matrices are defined by selecting L random 



vectors x q independently in the <i-dimensional hypercube [0, 1] . These vectors are called the vectors 
of phases. Given a scale parameter N, consider the function defined by 7 : {0, 1, . . . , N — l} d — > 
{0,1, . . . ,N d - 1} such that for every vector of integers m = (4,4, • • • 4d) G {0, 1, . . . , N - l} d 
the value 7(111) is equal to 

d 

7(m) : J>> % 
i=i 

It is easy to see that this function is a bijection over the set {0, l,...,N d - l}. Now we define the 
N d x L matrix as 

v !!H 9 ) := i ex p( 2fi ( m ' x 4 ( L2 ) 

For the case d = 1 we drop the upper index and we denote these matrices by V, for d > 2 we 
denote them by Y^ d \ 

Random Vandermonde matrices and their extended versions are a natural construction with a wide 
range of applications in fields as diverse as finance [12) . signal processing [13], wireless commu- 
nications [14j . statistical analysis j3], security [18] and biology [19] . This stems from the close 
relationship that unit magnitude complex Vandermonde matrices have with the discrete Fourier 
transform. Among these, there is an important recent application for signal reconstruction using 
noisy samples (see |13j ) where an asymptotic estimate is obtained for the mean squared error. In 
particular, and as was shown in |13j . generalized Vandermonde matrices play an important role 
in the minimum mean squared error estimation of vector fields, as might be measured in a sensor 
network. In such networks, the parameter d is the dimension of the field being measured, L is the 
number of sensors and iV can be taken as the approximate bandwidth of the measured signal per 
dimension. This asymptotic can be calculated as a random eigenvalue expectation, whose limit 
distribution depends on the signal dimension d. In the case d = 1 the limit is via random Vander- 
monde matrices. As d — > 00 the Marchenko-Pastur limit distribution is shown to apply. Further 
applications were treated in [T7] including source identification and wavelength estimation. 

One is typically interested in studying the behavior of these matrices as both and L go to 
infinity at a given ratio, limAr^oo = j3. In [T7], important results were obtained for the case 
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d = 1. In particular, the limit of the moments of V*V was proved and a combinatorial formula 
for the asymptotic moments was given under the hypothesis of continuous density. In [21], these 
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results were extended to more general densities and it was also proved that these moments arise as 
the moments of a probability measure supported on [0, oo). This measure depends of course 
on the measure v, the distribution of the phases, and on the constant j3. 

In [21], the behavior of the maximum eigenvalue was studied and tight upper and lower bounds 
for this one were found. Here we extend these results and study the maximum eigenvalue of the 
d-fold extended Vandermonde matrix. More specifically, we study the asymptotic behavior of the 
maximum eigenvalue of the matrix V^-* \~( d ) and derive upper and lower bounds. 

A natural question is to examine the behavior of the smallest singular value as N — > oo, and this 
paper is one of the first to address this question. Here we restrict to the case d = 1. On one hand, 
it is clear that if L > N then the matrix V*V is of size L x L and rank N. Therefore, if f3 > 1 the 
asymptotic limit measure has an atom at zero of size at least /3 — 1. On the other hand, if L = N 
then with probability one the random matrix V*V has determinant equal to 

det(V*V) = | det(V)| 2 = • 11 \e 2 ^ - e 2m ^\ 2 . (1.3) 

1<P<9<AT 

This determinant is zero if and only if there exist distinct p and q such that 9 P = 9 q . This is an 
event of zero probability if the probability measure has a density. Therefore, the minimum singular 
value Ai(iV) is positive with probability one and converges to zero as N increases. In this work, we 
show that with high probability Ai(iV) < N 2 exp(— C\fN). As a consequence of our argument we 
show that with high probability 

max | J[\z - Zi\ 2 : \z\ = lj > exp(Cv / iV), (1.4) 

where z k = e 2lTi6k and {6>i, . . . , 6 N } are i.i.d in [0, 1]. Moreover, we explicitly determine the constant 
C. We believe that this may prove to be of independent mathematical interest. Additionally, we 
show the absence of finite moments for the matrix (V*V) _1 . 

Lastly, we present a generalized version of the previously discussed random Vandermonde matrices. 
More specifically, consider an increasing sequence of integers {k p } p % 1 and let {6±, . . . , On} be i.i.d. 
random variables uniformly distributed in the unit interval [0,1]. Let V be the N x N random 
matrix defined as 

V{i,j):=^=zf (1.5) 



where Zj := e 5 ' j . Note that if we consider the sequence k p = p — 1 then the matrix V is the usual 



random Vandermonde matrix defined in Equation (1.1). We study the limit eigenvalue distribution 
for the matrices X := VV* and in particular their asymptotic moments. In this work we find a 
combinatorial formula for its moments and we show that for every sequence there exists a unique 
probability measure on [0, oo) with these moments. Finally, we show that for the sequence k p = 2 P 
the limit eigenvalue distribution is the famous Marchenko-Pastur distribution. 

The rest of the paper proceeds as follows. In Section [2j we present some preliminaries in random 
matrix theory, the Littlewood-Offord theory, we set up some notation and terminology and we 
review some known results for random Vandermonde matrices. In Section [3j we prove a formula 
for the trace log and log determinant of the random Vandemonde matrices. We also prove the 
absence of finite moments for the matrix M*M where M = V -1 . In Section |4j we show upper 
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and lower bounds for the behavior of the maximum singular value of V" in the general case. In 
Section [5j we study the behavior of the minimum singular value of V. In Section [6j we present 
some numerical results that suggest the absence of an atom at zero for the limit eigenvalue for the 
square case. In the last Section, we analyze the moments and limit eigenvalue distributions of the 
generalized version of the random Vandermonde matrices as described before. 



2. Preliminaries 



2.1. Random Matrix Theory Preliminaries. Throughout the paper we will denote by A* the 
complex conjugate transpose of the matrix A and the matrix Ijy will represent the N x N identity 
matrix. We let Tr be the non-normalized trace for square matrices, Tr(A) := 5^i=i an, where an 
are the diagonal elements of the matrix A. We also let tr^r be the normalized trace, defined by 
trjv(A) = -^Tr(A). Let A^v = (°ii( w ))iy=i be a random matrix where the entries ajj are random 
variables on some probability space. We say that the random matrices Ajv converge to a random 
variable A in distribution if the moments of A^v converge to the moments of the random variable 
A. We will denote this by A at — > A. 

One should notice that for a Hermitian N x N matrix A = A*, the collection of moments corre- 
sponds to a probability measure /Ua on the real line, determined by trAr(A fc ) = J R t k dfj<ji(t). This 
measure is given by the eigenvalue distribution of A, i.e. it puts mass on each of the eigenvalues 
of A (counted with multiplicity): 

1 N 

MA = ^£fci, (2-1) 
i=l 

where Ai, . . . , Ajv are the eigenvalues of A. In the same way, for a random matrix A, fix is given 
by the averaged eigenvalue distribution of A. Thus, moments of random matrices with respect to 
the averaged trace contain exactly the type of information in which one is usually interested when 
dealing with random matrices. 

Consider an N x L random Vandermonde matrix V with unit complex entries, as given in Equation 



(1.1). The variables 6{ will be called the phase distributions and let v be their probability distribu- 
tion. It was proved in [T7] that if dv = f{x) dx for f(x) continuous in [0, 1] then the matrices V* V 
have finite asymptotic moments. In other words, the limit 



m 



W = lim E 



iV->-oo 

exists for all r > 0. Moreover, 



tr L (V*V) 



(2.2) 



t (/3) 



where K P)ll are positive numbers indexed by the partitions p. We call these numbers Vandermonde 
expansion coefficients. 

The fact that all the moments exist is not enough to guarantee the existence of a limit probability 
measure having these moments. However, it was proved in [21] that the eigenvalues of V*V 
converge in distribution to a probability measure fj,p jU supported on [0, oo) where (3 = limTv-s-oo jj- 
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More precisely, 

= I fdfi^(t). 
Jo 

In [21], the class of functions for which the limit eigenvalue distribution exists was enlarged to 
include unbounded densities and lower bounds and upper bounds for the maximum eigenvalue were 
found. We suggest the interested reader to look at the articles [T7] and [2T] for more properties on 
the Vandermonde expansion coefficients as well as methods and formulas to compute them. 

2.2. Littlewood-Offord Theory Preliminaries. Let v\,...,v n be n vectors in which we 
normalise to all have length at least 1. For any given radius A > 0, we consider the small ball 
probability 

p{v x , . . . ,v n , A) := supP(r/iUi + . . . + rj n v n 6 B) 
B 

where rjx,...,rj n are i.i.d. Bernoulli signs (i.e. they take values +1 or —1 independently with a 
probability of 1/2 of each), and B ranges over all (closed) balls of radius A. The Littlewood-Offord 
problem is to compute the quantity 

p d (n,A):= sup p(vi, ...,v n ,A) 

where vi,...,v n range over all vectors in M rf of length at least one. Informally, this number measures 
the extent to which a random walk of length n (with all steps of size at least one) can concentrate 
into a ball of radius A. 

The one- dimensional case of this problem was answered by Erdos. First, one observes that one can 
normalise all the Vi to be at least +1 (as opposed to being at most —1). In the model case when 
A < 1, he proved that 



Pi(n,A)= I , " )/2» 



2 +o(l) 



in/2} J' 

when < A < 1 (the bound is attained in the extreme case v\ = . . . = v n = 1). A similar argument 
works for higher values of A, using Dilworth's Theorem instead of Sperner's Theorem, and gives 
the exact value _ 

p l( n,A) = W H W= V* ( 2 .4) 

j=l \ 3/ V 

whenever n> s and s — 1 < A < s for some natural number s, where ("),...,(") are the s largest 

— — ' \mi/' ' \m a / ° 

binomial coefficients of (?),... , (") . See the reference |20j for more details on the Littlewood-Offord 
Theory. 

3. Trace Logarithm Formula and the Inverse of a Vandermonde Matrix 

3.1. Inverse of a Vandermonde Matrix. Given a vector x in C N , we define cr™(x) to be the 
sum of all r-fold products of the components of x not involving the m-th coordinate. In other 
words, 

<=y: n 
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where p™ is a subset of {x\, X2, ■ ■ ■ , x m -i, x m +i, ■ ■ ■ , X]y} , with cardinality r. 

The following Theorem was proved in 

Theorem 1. Let V be a square N x N matrix given by 

V := 



X-2 



1 

X N 



\ 



V x 



N-l 
1 



JV-1 



.N-l 



with no entries. Then its inverse M := V is the matrix with entries 



M(m, n) 



(-1) 



N—n„m 



N-n 



X) 



n 



Xj ) 



with m, n G {1, 2, . . . , TV}. 



Remark 1. Let v\ < v-i < . . . < un be the eigenvalues o/M*M and let X± < A2 < 
corresponding eigenvalues o/V*V, which are the same as for VV*. Note that 

-1 



i/fe = A 



jV-(jfc-l) 



(3.1) 



< \m be the 



and in particular vn = A 1 . Therefore, to study Ai it is enough to understand the behavior of v^. 



Here we prove a Theorem about the trace log formula for the random Vandermonde matrices and 
a Theorem about the non-existence of the moments of M*M, but first we need a Lemma. 

Lemma 1. Let M be an invertible N x N matrix with columns X\, . . . , Xn and let R\, . . . , Rn be 
the rows o/M _1 . Let Vj 6e i/ie subspace generated by all the column vectors except Xi, i.e., 

Vi = sp&n{X 1} X 2 , . . .,Xi-i,X i+ i, . . .,X N }. 

Then the distance between the vector Xi and the subspace Vi is equal to 

dist(X h Vi) = tt^t. 

H-n-ill 

Moreover, 

N N 

^A^M^M)- 1 =J2dist(Xi,Vi)- 2 . 

8=1 8=1 



Proof. The result follows from an identity involving the singular values of M. By definition, the 
inner product {Rk,Xf) = 5k/ so that Rk is orthogonal to W Hence, 



d(X k ,V k ) = 




Let A fe (M*M) be the eigenvalues of M*M. Then, 

N 

Tr ((M- 1 )*!^- 1 ) =^A 4 (M*M)- 1 . 

i=i 
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On the other hand, 



N 



Tr (M-')'M- 1 = £ l(M- 1 ) M | a = EHB*ll 2 - 



Kkl<N 



k=l 



□ 

Theorem 2. Let V be a square random Vandermonde matrix of dimension N with i.i.d. phases 
distributed according to a measure v with continuous density f(x) over [0, 1]. Then 



E 



trjvlog(V*V) = (N- l)E(log|l -e 2m0 \) -log(iV). 



(3.2) 



Proof. The matrix V is an N x N random Vandermonde matrix as in Equation (1.1), with i.i.d. 
phases distributed according to a measure v with continuous density f(x). Let < Ai < A2 < 
. . . < Atv be the eigenvalues of V*V. Note that Ai > with probability one. It is clear that, 

JV 

logdet(V*V) = ^log(A i ). (3.3) 



i=l 



On the other hand, 



hence 



det(V*V) 



N N 



n 



\ e 2me p _ e 2me q |2 



l<p<q<N 

logdet(V*V) = ^21og(|e 2 ^ - e 2nie *\) - Nlog(N). 



(3.4) 



p<q 



Since the phases are identically distributed it is easy to see that the expectation E[logdet(V*V)] 
is equal to 

N(N — 1) 



E 



21og|e 2 ^ - e 2ni6 "\ 



Nlog(N) 



which is 



N(N — 1) E 



log 1 1 — e 



27ri0| 



iVlog(iV). 



(3.5) 



(3.6) 



Since for every invertible Hermitian matrix A we have that Trlog(A) = logdet(A), in particular 
we have that 

t rA rlog(V*V) = llogdet(V*V). (3.7) 



N 



Combining equations (3.5) and (3.7) we see that 



E 



tr w log(V*V) =(iV-l)E[log|l-e 2 H] -log(JV). 



(3.8) 
□ 



Theorem 3. Let V be a square N x N random Vandermonde matrix. Then for every N > 2 the 
matrix V* V is invertible with probability one and 



E(trjv((V*V)-P)) =00 



for every p > 1 . 
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Proof. It is enough to prove the case p = 1 since the other cases follow from this one. From Lemma 
[Hwe see that 

1 N 

E(tr iV ((V*V)- 1 )) = - £ E(|M(m,n)| 2 ) 



N 

m,n=l 



where M = V .In particular, 



E(tr^((V*V)- 1 ))>lE(|M(l,iV)|2). 
Now using Theorem [T] we see that 

JV-1 



E(|M(1,JV)| 



2\ 



1 f 2n d6 



27T./ \l-e l 



2 



OO. 



□ 



4. Maximum Eigenvalue 

Let X be the Lx L matrix defined as X := V^ 2 ) V^ 2 ). Suppose that the phases (6k,ipk) are selected 
i.i.d. on [0, l] 2 . We further suppose, as in [21], that (9k,ipk) — (9i, V'l) given 9x,ipi has a conditional 
density (which exists for all (9i,tpi) if it exists for one) and denote this density by / ignoring its 
dependence on (6*1,^1) as it will only appear through ||/||oo- It can be shown that X has the same 
eigenvalues as the matrix A whose entries are 

A(k, m) := D N (2Tr(9 k - 9 m )) ■ D N (W^* - V™)) (4.1) 

where 

sin(^x) 
d n{x) := . 

iVsm(§) 

is the Dirichlet kernel (see [21] for instance). Similarly, in the case d > 3 we obtain a product 
of d Dirichlet kernels. Subsequently, A will be used to construct upper and lower bounds for the 
maximum eigenvalue Xl. 

We now proceed to obtain asymptotic upper and lower bounds for d-fold Vandermonde matrices. 
We will first focus on the case d = 2 and retain the notation of Section [2j In what follows we will 
prove the following Theorem. 



Theorem 4. Let\( d > be the N d x L d-fold random Vandermonde matrix defined in (1.2). Let Xl 
be the maximum eigenvalue of the matrix (V^)*(V^). Then, if N, L — > 00 such that [3 = jh € 
(0, 00) there exist constants C e such that for all C > C e 

F(\ L > Clog(^) + „) < e-« ^ (4.2) 



where 5 = C — C e . 
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Proof. The line of argument follows that of [2T] . We begin with the following upper bound for the 
Dirichlet function proved in |21j . 



' N I 



Asin(M) 



N/2 



k=i L 



2ir(fc-l) 2irk 
N ' N 



(W). 



(4.3) 



To apply the bound, let p a fi be the probability with a,b £ Z such that 

6k — Oi £ [(a-l)/N,a/N\ 

and 

^-Vi e [(6-i)/a,6/at]. 

Define, 

Qa,b ■= Pa,b + P-a,b + Pa -b + P-a-b 

Then, it is easy to see that 



Q a ,b < 



C a 



(4.4) 



where the function / is the probability density over the square [0, l] 2 . Next for the magnitude of a 
term in the first row, corresponding to 9\ , ipi we find that 



N/2 N/2 



N/2 N/2 



where X is defined as 



a=l 6=1 



a=l 6=1 



t 

' ab 



(4.5) 



X 



D N (27r(9 k - ft)) • D N (2^ k - 4,!)) 



for k 1 and the upper bound does not depend on k or (6\ If R is any row sum of the entries 
of the matrix X it follows that 

1 L-l 



E[exp(i?)] < e 



1 + C 1~fp~ H N/2 



< e 



(3 N C q (e-l)H 2 N/2 



(4.6) 



by taking t = 1 and using the fact that (1 + x/A) < exp(x) and that /3at — )• /3 as iV — > oo where 
■Hp := Y%=i 1A- Since 



F Ar = logiV + 7+^ + 0(iV- 2 ) 



; A 



(4.7) 



it follows that 

E(exp(E)) < e ^(iogA) 2 =iY Ceio gJ 

for some suitable constant C e > 0. Applying the union bound to the maximum row sum Y and 
then, Markov's inequality with C > C e + 5, we get 



y > C(logA) 2 + u) < e~ u 



A 2 



, jyC e log A _ g — U 



A 2 



jyCTogA ° ^5 log A' ( 4 ' 8 ) 

Since the maximum eigenvalue is upper bounded by the maximum row sum of magnitudes, then 
y > A/,. This concludes the proof for the case d = 2. 

For the case d > 2 one obtains a (i-fold product of Dirichlet functions so that the exponent of the 
Harmonic function H is d instead of 2 in (4.6) and also in (4.4). Constant terms are also suitably 
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modified and with these changes carried through to (4.8) where 2 is again replaced with d, the 



result follows as before. □ 

The following Corollary is stated without proof. 
Corollary 1. There exists a positive constant B such that 

IE [Xl] < B \og{N d ) + o(l). 

4.1. Lower Bound. The purpose of this subsection is to present the following lower bound for 
the maximum eigenvalue A^, in which we suppose that phases are provided according to a joint 
continuous density / bounded away from 0. 

Theorem 5. Let Xl be the maximum eigenvalue then if N, L — > oo such that linijv_> 00 ^ = /3 € 
(0, oo) there exists a constant K such that 

/ loe N d \ 

F[X L >K - f a . ) =l-o(l). (4.9) 
V loglogiVV y J y ' 



Proof. We rely on the equivalent matrix given in (4.1), which is an L x L matrix. For notation 
simplicity we specialize to the case d = 2 since the case d > 3 follows a similar argument. To 
construct our lower bound by analogy with the 1-dimensional case, we want to obtain large numbers 
of points (6, ip) which lie close together, that is points such that \9k — m \ and \ipk — ip m \ are both 
small. By hypothesis the joint measure has a continuous density / bounded away from 0, throughout 
[0, l] 2 . Therefore, it follows that 

rj := mm {f (9,^) : (0,^) G [0, l] 2 } > 0. 

We randomly select L phase vectors in [0, l] 2 according to this measure. Let us divide the square 
[0, l] 2 in kN 2 equal squares. Each of of these squares will receive at least ^2?7 points on average. 
This is an occupancy model and we are interested in the square which has the maximum number 
of points. The number of such points is at least 

log log N 1 

with probability 1 — o(l), for all a € (0, 1) independently of the mean number of ball per urn (see 
for instance |15j). By dividing more finely and then taking N sufficiently large the entries of the 
matrix can be made arbitrarily close to 1 and so have maximum eigenvalue arbitrarily close to 
k(N 2 ). Finally, using the interlacing Theorem as done in [21] on A the result is obtained as a lower 
bound on the maximum eigenvalue. □ 



5. Minimum Eigenvalue 
In this section we focus on the behavior of the minimum eigenvalue Ai for the case d = 1. Consider 



the matrix A as in Equation (4.1) and all its 2 x 2 principal sub- minors. These matrices are 
symmetric and the minimum eigenvalue for the sub minor determined by phases Ok and 0£ will be 
denoted X^/- In other words, X^i is the smaller root of the equation 

(1 - A/,, r = D N (2ir(9 h - " ■ l2 
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Taking square roots and again applying the interlacing Theorem we obtain, 

Ai(iV) < min [l - D N (2ir(0 k - t )) 

Let 7 and Cat be defined as 7 := N 2 min^ \6 k — 9(\ and 

C N := min (l - D N (2n(9 k - e )j\) 

From a result of de Finetti (see [9] for a reference) , 

F(^mm\e k -0 £ \ >S^=(1-N6) 



N-l 



where := max{x,0} and this operation is taken first. Substituting rj/N 2 and taking the limit 
as N — > 00 we obtain that 

V 



lim P min \0 k - 9A > 

N^oo \ k,£ 1 N 2 

On taking the Taylor expansion of to second order around the origin we obtain 

2 2 

min (l - D N (2n(9 k - 9 e )j) = ^ + o(N- 2 ). 
Therefore, the following limit holds, 

2 

lim F(C N < ) = 1 - e" 7 . 

Then, it follows that Ai(iV) < 0(iV~ 2 ) as iV -> 00. As we will now show, the approach to zero is 
much more rapid. To obtain better estimates for \i(N) as — > 00, we now consider the maximum 
eigenvalue of the inverse matrix. 

Given an N x N matrix M, we have the following matrix norms ||M||i := sup^- (X^a-=i l-^fcjl) an d 
||M||2 := \J Aj\r(M*M). The following inequality is well known (see [10] for more details), 

— L||M||i < ||M|| 2 < ||M||i. 

We now prove the following Lemma which will be used later. 

Lemma 2. Let A(z) = aQZ n +aiz n ~ 1 +. . .+a n be a complex polynomial and let A* := max| 2 | =1 |^4(z)| 
be its maximum in the unit circle. Then 

\a n \ + . . . + \a \ .* i . 

— <A < \a n \ + ■ ■ ■ + \a \. (5.1) 

n + 1 

Proof. The second inequality follows immediately from the triangle inequality, so we concentrate 
on the first. It will be enough to show that 

\a k \ < A* (5.2) 
for all k. By applying Cauchy's integral Theorem and by using the fact that 

z~* dz = 

|*|=1 
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for alH ^ 1, it is obtained that, 
for all k. Therefore, 



1 f A(z) , , , 



2?r Jo 



where the first inequality follows by upper bounding the integral, taken as a line integral with 
respect to 9 around the unit circle, and the second inequality follows from the definition of A*. By 
applying the inequality to each k we obtain the required lower bound. □ 

In the following steps we find a bound on Ai(iV) in terms of the maximum of a polynomial with 
roots on the unit circle. We begin with some definitions. Let Zk = e 2m9k be the values determining 



the random Vandermonde matrix as in Equation (1.1). Let P(z) be the polynomial defined as 

N 

P{z) := J] (z - z k ) . 

k=l 

We further denote 

| Z Zp 

Let M = V- 1 be the inverse of the random Vandermonde matrix and let M(p, q) denote its entries. 
Define 

N 

P p :=J2\M(p,q)\. (5.4) 

q=l 

By Theorem [TJ we know that 

h ■= Trnr-n ( k ° 1 + ■ ■ ■ + |(T M • (5 - 5) 

llq^p \ z p ~ z q\ v 7 

In addition, let 



It follows from (5.5) and (5.1) that 



< max 
N ~ \z\- 



z(yN-\T p (z)\) <P P 

for all p = 1, . . . , N. The following Lemma is a direct consequence of the Hadamard's inequality. 



Lemma 3. Let zi,...,Zn be distinct points on the unit complex circle. Then there exists po 6 
{1, . . . , N} such that, 

U\z P0 -z q \<N. (5.7) 

We are now in a position to prove the following Lemma, which provides upper and lower bounds 
on the minimum eigenvalue of a Vandermonde matrix, in terms of the polynomial T p , defined in 



(5.6). 
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Lemma 4. Let Xi(N) be the minimum eigenvalue of the random Vandermonde matrix as defined 
in (1.1) and T p {z) be as defined above. Then, 

< A i(^V) < p 1 i^, mot ( 5 ' 8 ) 



N 3 maxp {max| z | =1 |T p (z)| 2 } max p {max| 2 | =1 |Tp(z)| 2 } 

Moreover, 



N 



2 



Xt(N) < = (5.9) 

max| 2 | =i n^ po k- " '- 



\z - Z q \ 



Proof. Since | |M| |i = max {(3 p : p = 1, . . . , N}, it follows that 

IIMI 



N 

On the other hand, 



- < VN ■ max I max \TJz)\ \ < \ \M\U. 

V \\z\=\ "J 



1 N 

rjTT < Xl(N) < 

M 2 ~ v 7 ~ M 12 



U 

from which we can deduce that 

m 7 1 \rr< M2\ ^ MX) < J 1 |Tn|2 v (5-10) 

JM^ maxp |maxi z | =1 max p |max| 2 | =1 |i p (z)| z j 

Using Lemma [3l we know that there exists po such that Ilg^pokpo ~ z q\ < N. We thus obtain that 



max \T po (z)\ > max ( — TT \z - z q \ | . (5.11) 

Z=l 2=1 I iV J - J - 



Therefore, using the fact that 
and Lemma 0] we see that 



max max |T p (z)| > max |T Po (z)| 

P \z\=l \z\=l 

N 2 

X 1 (N) < V2 . (5.12) 

max k | =1 \[ q ^ po \z - z q \ 2 

□ 



5.1. Stochastic Construction. Before stating our upper bound for the minimum eigenvalue, we 
introduce the following definitions. First, we define a random sequence via a realization of the 
Brownian bridge W° on [0, 27r], which satisfies VF°(0) = W°(2ir) = 0. A <p shift of the Brownian 
bridge is defined by 



W°{0) :-- 



W°((p + 9) - W°(ip) if < < 2vr - 

W°((p + 9 -2tt) -W°(ip) if 2vr - 99 < 6 < 2vr. 



Further, define the infinite sequence $ := {(f r }r>o to be the sequence of dyadic phases on [0,2-7r]. 
Given a realization of the Brownian bridge define the following function, 

r>27r 



Jo v 1 cos P 
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Figure 1. Realization of the Brownian Bridge and a Shifted Version. 











Aw 




\ ft 























Figure 2. J« for the Brownian bridge previous realization. 



for ip 6 [0, 27r]. Note that it is not clear that the above integral is well defined, since it may not exist 
behaves like 9~ l near and near 2tt. We will address this matter shortly. In 



as the fraction 



sin 9 
1— cos ( 



Figure [TJ we show a realization of the Brownian bridge and a shifted version with ip equal to 2>ir/2. 
In Figure [2j we show I v for the same realization as before. 

Using the sequence $ and the same realization of W° we construct the sequence of random variables 
I = {l r }r>0 as 

Ir ■ — I<p T ■ 

It will be shown that J„ is continuous on the interval [0, 2ir], and so there exists a value p* which 
determines the maximum value of J„, which we denote as I*. Since $ is dense on the unit circle, 
it follows that 

I* := sup {I r : r G N} (5.13) 

and its distribution is determined via the infinite sequence I r . We now show that the above random 
function is well defined and that integrals are a.s. finite. 

Lemma 5. Given a realization of the Brownian bridge W° , then a.s. the following integrals exist 
for all (p £ [0, 2ir) 



2t sin^ 



W°— 



cos ip 



dtp 



< oo 
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Proof. When <p> = we write the above integral as /. The Levy global modulus of continuity tells 
us that for standard Brownian motion B on [0, 2tt) satisfies almost surely 

r r \B(t + 8)-B(t)\ 
nm limsup —■ = 1 

s ±0 0<t<2n-8 w(d) 



where w(S) = y 28 log j (see [E] for a proof of this fact). Since W° is by definition, 

our argument will be the same no matter which value of tp is chosen because the Levy modulus 
applies to the entire sample path. We therefore set ip = 0. By definition of the Levy modulus, 
almost surely there exists 82 > such that 

|J?(t + *)-B(t)| 



w(S) 

for all < 8 < 82- Therefore, 



a(6) := \W°{iP + 8)- W°(tp)\ < 2w(5) + ^^-8. (5.14) 



2tt 



We may therefore split the integral as, 



2tt-5 2 sin?/; f' S ' 2 

w°w- — ^-#+/ w°w- 

1 - cos ip Jo 1 



sin?/; f 2 simp 



g 2 1 — cos ip J 1 — cos ip 

271 sinip 



dip 

Jo 1 ~ cus W 

J2tt~5 2 1- COS Ip 

The first integral is finite being the integral of a continuous function over the interval [82 , 2tt — 82} ■ 
We may further suppose that 62 has been chosen so that I^ tz^wH — 4 for < ip < 62 with the 
corresponding inequality in a similar neighbourhood of 2n. By choice of 82 we obtain that 







1 — cos ip 



"0 
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for sufficiently small 62- The same argument applies to the last integral. Since 10(62) gives a uniform 
bound the result holds for all ip € [0, 2tt). Continuity in ip follows by a similar argument, 

(•27T-.5 



lip I(p I 



< 
+ 
+ 
+ 
+ 



l-A-K — O 

J s (K - w $) r 



siwip 



\w°M)\ 



simp 



/ \w;m 

J2ir-6 

/Vkvoi • 

Jo 

/•27T 
J27T-5 



1 — cos ip 
sin^ 



cosV' 



(5.15) 



1 — cos ip 
simp 



dip 



cos ip 
sin ip 



dip 



1 — cos ip 



dip. 



Provided that < S < 62, the tail integrals are all at most 0(5 1 ^ 3 ) as before. We bound the first 
integral by two positive integrals, to obtain 



I-27T-S 

/ {w; - w$) 

J s 



simp 



-dip 



1 — cos ip 
< 2- sup \W°^)-W$(iP)\- 



simp 



dip 



0<</><2tt ' J 8 (1-COSVO 

< 2 • sup|^°(V>) - W$ty)\ • [(1 - cos iP)Y s 

< 12a(<5) (log 2 - log(l - cos 5)) 

provided \<p — <p\ < 5. Finally, 12a(<5) (log 2 — log(l — cos 5)) — > as 5 — > 0, which implies continuity. 

□ 



It therefore follows that I* is well defined. Let Tw(<p) be defined as 

1 1 N 

T N (<p) ■= -7= log \P(e^)\ 2 = — ^ log [2(1 - cos( 9 9 - 6 q ))] (5.16) 

V V q=l 

where P(z) is the random polynomial in the unit circle with roots {e t9q } q=1 as before. Furthermore, 
define the infinite sequence of random variables T^(ip r ) by evaluating the previous expression at the 
phases of <3>. Note that T/v cannot be defined as a random function in either C[0, 2ir] or in D[0, 2tt] 
as its discontinuities are not of the first kind. We remark that since there are only countably many 
<p r and the phases 6 q are i.i.d. and uniformly distributed then no <p r coincides with any 9 q almost 
surely so that the sum exists. We further observe that since, 

r2w 



and 

r 2ir 



r2ir 

/ log [2(1 - cos(V>))] dip = 
Jo 

f-2n 

/ log 2 [2(1 - cos(YO)] dip < OO. 
Jo 
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Note that the sequence Tjv(<^ r ) satisfies the central limit Theorem as a function of N for every 
r £ Z+. We consider the sequence TV := {Tjy(ip r )} r >o in the sequence space M°° with metric 



OG 



Po x, y) = > — ,2 

*—'l + \xe- vf\ 



e=0 ~ + \ x t-ye\ 

and using the ordering stated earlier. It is well known that this forms a Polish space ([!])■ We now 
derive one more Lemma for use later on. 

Lemma 6. Let Y be a function in D[0, 2tt] . Then Y is Lebesgue measurable, and its integral exists, 

Y(s)ds < oo. (5.17) 

o 

Furthermore, let Y n be a sequence of functions in D[0, 2ir] such that Y n — )• Y in D (i.e. with respect 
to the Skorohod topology) then 

"2tt r-2n 

Y n (s)ds^ / Y(s)ds. (5.18) 
J o 

Proof. The existence of the integral follows from Lemma 1, page 110 of |3J and the subsequent 
discussion which shows that functions in D on a closed bounded interval are both Lebesgue mea- 
surable and bounded. The former follows from the fact that they may be uniformly approximated 
by simple functions, a direct consequence of Lemma 1 and the latter also. 

Convergence follows from the Lebesgue dominated convergence theorem. This holds since the 
sequence Y n is uniformly bounded, by a constant so the sequence is dominated. Second Y is 
continuous a.e. with pointwise convergence holding at points of continuity, as a consequence of 
convergence in D see [3]. □ 

We now proceed to prove the following Theorem. 

Theorem 6. With the topology induced by the previous metric in M°° we have that the sequence 
T^v converges in distribution to the sequence I 

Tat^-I (5.19) 

as N — ^ oo. 

Proof. In order to do so we will Theorem 4.2 of |4j. Suppose that there is a metric space S with 
metric po and sequences Tat i£ , I e and Tjv all lying in S such that the following conditions hold, 

Tiv, £ Is (5.20) 

I £ I 

together with the further condition that given arbitrary rj > 0, 



limlimsupP(po(Tjv, e ,Tjv) > 7/J =0. (5.21) 

Then it holds that Tjv =>■ I. First, we define I e using a realisation of the Brownian Bridge as 
follows, 

I r , £ := [^Wlog2(l-cos^)]f- £ -^"V° rr ^^#. 
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The definition of the other sequence is more involved and so we defer it for a moment. We have 
shown that the limit integrals exist a.s. and so we only need to show that the first term converge 
to 0. Since log [2 (1 — costp) ] = O(loge) when e is small and in a neighbourhood of and 2ir we 
may invoke the Levy modulus of continuity, wrapped around at 2ir to obtain that this term is 

0(loge x o(e)) -»■ 0. 

Hence coordinate convergence of the integrals holds so that 

sin-0 



Ir.e + 



2tt- 



w° 



dip 







and it follows that I F 



1 — cos tp 

I as £ | 0, since the sign of the integral is immaterial. We have thus 



demonstrated the second condition of (5.21). Finally, we proceed by rewriting Tjv(y? P ) i n terms of 
the empirical distribution function F/v : [0, 2ir] — > [0, 1] determined by 



#{9 q :0<e q <^} 



N 



It is a straight forward observation to see that 

T N (cp, 



N log [2(1 - cos(y> r - i(,))]dF N (i/>) 







2tt 



log [2(1 



COS 'i 



\dF N>Vr $) 



where the change of variables, ip = ip — (p has been made. For ip G [0, 2ir) we define Fpf<p(ip) as the 
"cycled" empirical distribution function of Fn by 



#W<o q «p+i>} 



N 



FnA 2 ^) + 



#{0<6> g <V>-27r+y} 
N 



if <p < tp < 2tt — <p, 
if 2ir - ip < ip < 2ir. 



To define the sequence T/v, e (</iv) we split the integral into two parts as in £ and J Q e + and 
then use integration by parts on the first part, which yields the expression, 



TN,e{<Pr) 



N x 



27T-£ N 



2n-e / ^ 



sinip 



2tt J (1 — cos ip) 



dip. 



For later use we make the definition, 



W N „ := VN (f n ^) 



2vr 



This is not quite equal to the original sum, since 

r'2ir 



/ log [2(1 - cosip)]dip = 
Jo 



(5.22) 
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so that the ip terms do not give but rather cancel with p, £ to be defined in a moment. The 
remainder we express it as a sum, noting that we must include the mean, which is by symmetry, 



Me 



2 

2^ 



log [2 (1 — cos ip) ]dip = — (eloge — e + o(e)) . 



7T 



Define S s (tp) := {9 q : 9 q G [tp — e, ip + e]} and hence the sum may be written as 



1 



iv ^ 



log [2(1 - cos(w - 0,))] - VNn e . 



(5.23) 
(5.24) 



Denote the corresponding sequence as Zat j£ . Taking expectations we thus find that 

N 



E[Z Nje (<p r )] 



log [2(1 - cos(^))] # - Viv> 







is a sequence of random variables with mean. We finally write, 

Tat = Tjv.e + Zjv.e- 



We now proceed to demonstrate the first condition of (5.21), namely that, Ttv, £ 



(5.25) 
Since 



TjV,e(¥v) is a functional of an empirical distribution and therefore of a process lying in L>[0,2-7r]. 
Define the random sequence J £ defined for / £ D[0, 2tt], /(0) = f(2ir) = with the component 
term, 

sin^ - ,i2tt-£ 



It is well known that 



2-rr-e 



cos 



/^(V>)log[2(l-cos^)] 



(5.26) 



in D, which implies that 

w NiVr 

as N — > oo for all r. The result will follow on showing that J £ defines a measurable mapping 
J £ : D[0, 2vr] -> IR°° in L>[0, 2tt]. Since 

Js,r(W N ) = T N>e ((p r ), 

we may therefore apply Theorem 5.1, Corollary 1 of ^ which states that if Wjv =^ then 
•^e(WOv) Je(W°), (and hence Tat i£ =4> I £ ) provided that we verify 



W° € Dj 



0. 



(5.27) 



To deal with measurability question we first observed that the coordinate maps are measurable and 



61 that J e ^ r is measurable 
the sequence of integrals 



since sin^/(l — cos-i/O is continuous in [e, 2ir — e], it follows by Lemma 
for each r and hence so is the sequence mapping J £ . Again by Lemma [6 
convergence with respect to po- This leaves only the final term. However since the limit W° is 
almost surely continuous it follows almost surely that 

UQjr-e) -> W° t {2it-e) 

for each r if / — > W° in L>[0,27r]. Thus the corresponding sequence converges with respect to po 
also and so (5.27) holds. The proof of the first condition is concluded. 



It remains to demonstrate (5.21). Here we will use the union bound and Chebyshev's inequality. 
This is because the various Z^^ £ {ip r ) in the sequences are dependent, as they are determined via 
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the same 9 q . Nevertheless they are of course themselves the sum of i.i.d. random variables. In 
determining the variance, we may work with ip r = without loss of generality. The variance of one 



of the i.i.d. summands in (5.24) is determined as 



2tt 



log 2 [2(l-costy))]#-M6 



< oo. 



(5.28) 



Since for small e > we have log [2(1 — cos(VO)] = 0(2 log -0) + o(ip) the integral is a £ = 0(elog 2 e) 
as the integral of log 2 x is xlog 2 x — 2xlogx + 2x. It follows that cr 2 — > as e — > 0, which is the 
variance of the entire sum by independence and as it has been scaled. 



Now fix rj > 0. By definition of po and from (5.25) we obtain that, 



^ 1 + \Z N ,eWr)\ 

Let be such that YlrLii +i 2~ r < Now we apply the union bound to the remaining R v + 1 
summands to obtain that 



R 



\ZN,e{Vr)\ 2 - 

1 + 1^(^)1 



> r?/2) < £>(|Z JV , e ( ¥ v)|2- r > 



r=0 

< ^^(i^ + l) 2 



2 (^ + 1) 



(5.29) 



r=0 



r? 2 2 2r 



Hence, 



< a. 



limsupP A)(Tjv, e ,Tjv) >r/ < 

iV 



216(^ + 1) 



3r/ 2 

16(i^ + 1) V 



2^2 



37/2 



and the RHS goes to as e to 0, for each r\ > 0. Hence we obtain (5.21) as required. We have 
therefore verified all conditions and Theorem [6] is proved. □ 



We are now in a position to state our main result for an upper bound on the minimum eigenvalue 
of a random Vandermonde matrix. 

Theorem 7. Let \\(N) be the minimum eigenvalue of the square N x N matrix V defined in 



1.1). We further assume that the phases 9\, ... ,6 n are i.i.d. and drawn accordingly to the uniform 



distribution. Then 

Ai (iV) < 2N 2 exp ( - VNT^ /2) (5.30) 
where Tjy- := limsup r Tjv(^v). Moreover, given a > 0, 



limmfp(Ai(AT) < 2N 2 exp ( - \/iVa/2)) >f(i* >aj. (5.31) 
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Proof. Using the definition of T/v(</v) in (5.16) and the upper bound on the minimum eigenvalues 
as given in (5.9) as well as the inequality 



max 



k TT \z — z a \ > - maxTT 



£ — 2 n 



<i¥=po 



we obtain that, 



Ai(iV) < 2iY 2 exp(-v / iVr 7 v(^r)/2). 



Since this equation holds over all r E Z + it follows that, 

Ai(iV) < 2iV 2 exp ( - v / iVlimsupT A r(99 r )/2). 



(5.32) 



Since this holds for all N G N we obtain ( 5.30| ), which is the first part of the Theorem. Now define 
the random variable, 

2 X 1 (N) 
log 



and further for any given R define Tn,r '■= max r <^j {T/v(c^, r )}, and similarly define Then for 
any fixed R and a > 

^>a) >Pf3V, Ji >a^ 



by definition of Tn,r- By weak convergence, since the set is open, and by Theorem 2.1 of [I] we 
obtain, 

liminfpfLAT > a) > lim inf P f TV fl > a) > f(i r > a 



N 



N 



as a consequence of Theorem [6j Finally by almost sure continuity it holds that Ijj f I* almost 
surely, and so by the monotone convergence theorem we see that P^J* > a^j = lim^P^/|j > a^j 
which implies our result. □ 



5.2. Analytical and Combinatorial Construction. In this Section we present an analytical 
and elementary argument for the upper bound of the minimum eigenvalue. Let zi,z%, . . . ,zn be 
complex numbers in the unit circle and let P(z) = YiiLi ( z ~ z i) be the polynomial with these 
roots. We want to estimate max| z | =1 |-P(.z)| when the roots {zi}f =l are uniformly distributed i.i.d. 
random variable in the unit circle. 

Lemma 7. Given P{z) as before there exists \ w\ = 1 such that \P{w)P{—w)\ = 1. 

Proof. Consider the function ^(z) = log|P(z)| + log |P(— z)\. This function is continuous except 
at the values {zi,— z\, . . . , zjy, — zn} where it has a vertical asymptote going to — oo. Therefore, we 
can consider this function as a continuous function from the unit circle to [— oo, oo) with the usual 
topology. On the other hand, it is clear that Jj z | =1 ^( z ) = 0. Therefore, there exist w such that 
= and hence \P(w)\\P(-w)\ = 1. □ 

Consider the following construction. We first randomly choose the points {zi}fL 1 and we consider 
the set of pairs V := {(^i, —z\), . . . , (zjy, — zn)}- Note that by changing Zi by — Zi does not affect 
the value of the point w in the previous Lemma. Hence the set V determines the point w. Now we 
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fix this point and consider a, := \w — z%\ and ft := \w + Zi\. Since \P(w)P(—w)\ = 1 we see that 



ni=i a iPi = 1- It i s a l so clear that ft = -i/4 — af. 



Let y be the random variable defined as y : {1,-1} 

N 

y(v 1 ,v 2 , ...,v N ) = ^Vilogipti/Pi) 
t=i 

taking signs i.i.d. at random with probability 1/2. It is not difficult to see that E(y) = where the 
average is taken over the set {1, — 1}^. Note that y(v\, . . . , vn) is equal to 

N 

y(vi,...,v N ) = ^^(log(ai) - log(ft)) 
i=l 

= log | ^...^O) I - log|P(„ li ... i „ Jv) (-«;)| 
where P( Vl VN )(z) is the polynomial with roots ViZi 

N 

P(vi,...,v N ) (z) = n {z - ViZi) (5.33) 
i=i 

and w is as in Lemma [7J Since \P(vi,...,v N )(~ w )\ = \P(v 1 ,...,v N )( w )\~ 1 we see that 

y(vi,...,v N ) = log|P ( „ lt ... )t)JV .)(w)| 2 . (5.34) 

Theorem 8. Let 7 := log ( / /f\ ) • For every e > the following holds 



sin(7r/8) 

!/(«!,..., vjv) I >7v / ^ev / iV/2) > 1 - e. 

Proof. By changing 2$ to — z%, if necessary, we can assume without loss of generality that on > ft. 
The point u; is equal to ro = e ie for some phase in [0, 2ir). Let .A be the set 

A:={e** : G [0 + tt/4, + 3vr/4] U [0 - 3vr/4, - vr/4]}. 

The total length of the set A is n and hence the probability of random point z% to belong to A is 
equal to 1/2. On the other hand, it is easy to see that if z% belongs to A c and since by assumption 
on > ft we see that 

log(Vft) > log ( C -^r^r) =: 7 - 0.8814. 
^sin(7r/8) J 

Let us order the values of log(aj/ft) in increasing order. Up to a re-numeration we see that 
< . . . < log(a m //3 m ) < 7 < log(a m+ i//Wl) < • ■ • < log(a N /f3 N ). 

The value of m is a random variable that converges almost surely to N/2 as N — > 00. Without loss of 
generality and for notation simplicity we will assume that m = N/2, however, and as the argument 
shows this is not strictly necessary. Let e > and let us consider the value sup^- P(y(«i, . . . , vn) € I) 
where I ranges over all the closed intervals of length r y^/7^e^/N in the real line. By applying the 
Littlewood-Offord Theorem discussed in the preliminaries section we see that 

sup ^f(y(v) € I I conditioning on Vi for i < [N/2\^j j = e + o(N~ 1/2 ) (5.35) 
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for N sufficiently large. On the other hand, 

sup{p(y(i;i,...,i;jv) El)} 

is equal to 

2jv/2 SU P {p(j/( w ) E I I cond. on the first Vij } < e 

— ,«LjV/2j) * 



where the last inequality is coming from Equation (5.35). In particular taking I to be the interval 
I = [_ 7 ^ e ^/]v /2, j-^/neVN /2] we conclude that 

p(|y(i>i,...,^)| > 7V WiV/2) > 1-e. 

□ 

Theorem 9. Given e > we have that 



for N sufficiently large. 



max|P(z)| 2 > exp(7v / vrev / iV/2) ) >l-e (5.36) 

1 2=1=1 



Proof. As we mentioned before we start by randomly generating n pairs of diametrically opposite 
points (zi,—Zi), and then find w as in Lemma [7J Finally fix the Z{ by the N independent coin 
flips (ui, . . . , vjv). Now we condition on this event. We observed before that \P(vi,...,vn)(~ w )\ = 
\ p (v 1 ,...,v N )(w)\~ 1 - Therefore, 

log\P {vii ... !VN )(w)\ = -log|P ( „ 1) ... iUjv) (-u;)|. 

Let a[w) be 

a{w) = \og\P^ tVN) {w)\ -log\P( Vlr .. jVN )(-w)\ 
= i °s\P(v 1 ,...,v N )(w)\ 2 . 



Then by the previous Theorem and Equation (5.34) we see that 

\a(w)\ > 7^6^/2) > 1 - e. 
Since o(— w) = —a{w) we clearly see that 

P( max (a(w), a(-w)) > ^e^N /2) > 1 - e. 

Therefore, 

p(maxmax{log|P ( , li ... > , Jv) (z)| 2 ,log|P (l , li ... i , )Jv) (-z)| 2 } > > 1 - e. (5.37) 

Now removing the conditioning on the pairs {21, — zi, . . . , zat, —zn} we see that, 

P [ max|P(z)| 2 > exp {Jirey/N/2) J >l-e (5.38) 

V w=1 / 

for N sufficiently large. □ 
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Proportion of e-values smaller than 1 0" p for N=1 000 and 50 trials 




Value of p 



10 11 12 13 



Figure 3. Graph of the average proportion of eigenvalues smaller than 10 p as a 
function of p for N = 1000. 



Since we already saw in Lemma [4] that 

Al ^max„ |=1 |P(*)| 2 
the following Theorem follows immediately. 
Theorem 10. Given e > we have that 

'Ai < A^ 2 exp(-^re\/iV/2)) > 1 - e. 



(5.39) 



6. Numerical Results 

6.1. Atom at 0. In this subsection we present some numerical results for the behavior near the 
origin of the limit probability distribution of V*V. Let V be a square N x N random Vandermonde 
matrix with phases 9\, 62, ■ ■ ■ , On which are i.i.d. uniformly distributed on [0, 1]. We know that the 
empirical eigenvalue distribution of V*V converges as N — > 00 to a probability measure fi. One 
question that we would like to address is: does the measure \x have an atom at zero? 

Let {\i\f =x be the eigenvalues of V*V. Given e > let us denote by Gjq{e) the average number of 
eigenvalues less than or equal to e, i.e., 



G N (e) :-- 



N 



-E 



{Ai < e : i = 1,...,N} 

Therefore, if there is an atom at zero for the measure \x with mass n{0} = f3, the following holds 

(3. (6.1) 



inf liminf Giv(e) 

e>0 JV-s-oo 



In Figure [5J we plot Gn(W~ p ) as a function of p for N = 1000. These graphs suggest that if there 
is an atom its mass has to be smaller than 0.06. Further simulations suggest the absence of an 
atom at zero. However, at the moment, we are unable to prove this result. 
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Figure 4. T N (ip) for 2 15 dyadic rationals with N = 10 6 . 



0.7p 
0.6- 




Figure 5. max r Tjv^yy) for 2 15 dyadic rationals with N = 32 and 100,000 trials. 

6.2. Results for the Minimum Eigenvalue. Our theoretical results show that the minimum 
eigenvalue is decreasing extremely fast. In Figure |4j we show a realization of T/v for a subset of the 
dyadic fractions and a large value of N, which indicates its qualitative behaviour. The maximum 
value over the values of ip selected is around 4. In Figure [5j we have obtained an empirical pdf for 
the above maximum. Finally, we show the results for max r Tn n((p r ) when R = 2 15 . 

Finally, we present some numerical results for the behavior of the maximum of a random poly- 
nomial in the unit circle in the context of Theorem [9} In Figure [6j we show the graphs of 
21ogmax^ =1 | |P(#)| and ^fne^/N /2 as a function of N. This graph suggests that Theorem [9] 
could be slightly improved. 

7. Generalized Random Vandermonde Matrix 

In this Section we present a generalized version of the previously discussed random Vandermonde 
matrices. More specifically, consider an increasing sequence of integers {k p }JL 1 and let {#1, . . . , 9^} 
be i.i.d. random variables uniformly distributed in the unit interval [0,1]. Let V be the N x N 
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Similations vs Lower Bound 




20 40 60 80 100 120 140 160 180 200 



Figure 6. Graphs of the average of 2 log max^ =1 | \P(z)\ (blue) and 
(red) as a function of N where 2 logmax| z=1 | |-P(z)| was averaged over 1000 realiza- 
tions. 



random matrix defined as 



where zj := e 2w K Note that if we consider the sequence k p = p — 1 then the matrix V is the 



V{i,j):=-±=z* (7.1) 
ViV 



usual random Vandermonde matrix defined in Equation (1.1). We are interested in understanding 
the limit eigenvalue distribution for the matrices X := VV* and in particular their asymptotic 
moments. Let r > and let us define the r-th asymptotic moment as 



m r := lim E 

N^oo 



tr N (X r ) . (7.2) 



These moments, as well as the limit eigenvalue distribution, depend on the sequence {k p } p= i- 

Remark 2. It is a straight forward calculation to see that m,Q = mi = 1, m.2 = 2 and = 5 
no matter what is the sequence {k p } p *L 1 . The first interesting case happens when r is equal to 4- 
These is because r = 4 is the first positive integer where there is a non-crossing partition, namely 
the partition p = {{1, 3}, {2, 4}}. 

The next Theorem shows a combinatorial expression for the moments as well as the existence of 
the limit eigenvalue distribution. 

Theorem 11. Let {k p } p *L l be an increasing sequence of positive integers. Then 

m r = K P ( 7 - 3 ) 
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where V(r) is the set of partitions of the set {1,2, ... ,r} and 

K := Urn |S ff| , 
where \p\ is the number of blocks of p and 

S PiN :={{p 1 ,...,p r )e{l,2,...,NY : E k Vi = E k n+1 ), 



(7.4) 



(7.5) 



where Bj are the blocks of p. Moreover, there exists a unique probability measure p supported in 
[0, oo) with these moments. 



Proof. Given r > then 

tl-Ar(X r ) = — X {VUP2)X{p 2 ,Pd)---X{jp r ,pi) 



(pi,-,Pr) 



Tr+l 2—/ 2-^ h ' '• 



y (kp\ kp 2 ) (kp2 kp^) „(.kpr kpl) 

(pi,...,p r ) («l,...,ir) 



The sequence (ii,i2, ■ ■ ■ ,i r ) 6 {1> 2, . . . , N} r uniquely defines a partition p of the set {1, 2, . . . , r} 
(we denote this by {i\, . . . , i r ) h-» p) where each block Bj consists of the positions which are equal, 
i.e., 

Bj = {w jl ,...,Wj ] } 



where i,, 



. Denote this common value by Wj. Then 



IpI 



Tr+l e En 



Eigs, ( k Pi k p i+1 ) 



^ ^ -I- -I- 

(n,...,j r )i-s>p (pi,...,p r ) fc=i 



(7.6) 



Taking expectation on both sides we observe that 



E 



Eigs fc ( k Pi k Pi+i) 



if and only if YlieB k (^p» ~~ ^Pi+i ) = 0- Let S^jv be the r-tuples (pi , . . . , p r ) which solve the equations 

E k n = E W (7-7) 

i&B k ieB k 

for all the the blocks k = 1, 2, . . . , \p\ and let IS^jvI be its cardinality. Let K p be defined as 

K = lim l 5p,jy l 



The it follows from Equation ( 7.6 ) that 



m r = E ^p- 



It is straight forward to see that the set of solutions of (7.7) has r + 1 — \p\ free variables since one 
of the equations is redundant (the sum of all the equations sums 0). Therefore, for every partition 
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p the value of K p satisfies < K p < 1. Then the moments are bounded by the Bell numbers 
B r = \V{r)\. Define, 



— — r 
(3 r := inf m^ fc < B^ k < inf yk = \Jr. 

k>r k>r 



Hence, (3 r 1 > r l / 2 and therefore 

+oo 
r=l 

Therefore, by Carleman's Theorem [2] there exists a unique probability measure p supported on 
[0, +oo) such that 



roc 



m r = t n dp(t). 
Jo 

In other words, the sequence m r is distribution determining. □ 
Proposition 1. Let p € V(r) then K p = 1 if and only if the partition p is non-crossing. 

The proof of this results follows similarly to the one presented in [T7] for the sequence k p = p — 1 
and we leave it as an exercise for the reader. 

Example 1. Let r = 4 and let p = {{1, 3}, {2, 4}}. Then 

K - lim ^1 

where 

S p ,n = {{pi,P2,P3,Pa) e {1,2, . . . , iV} 4 : k pi + k P3 = k V2 + k P4 }. 

For the case k p = p — 1 it was observed in [T7] that K p = 2/3. As a matter of fact, it is not difficult 
to see that K p is the volume of the polytope 

K p = vol({(x, y,z) e [0,1] 3 : 0<x + y-z< 1}). 

This polytope is shown in Figure^ For the case k p = 2 P we see that 

S p ,n = {(pi,P2,P3,P4) € {1, 2, . . . ,iV} 4 : 2" 1 + 2^ = 2^ +2^}. 

For positive integers {a, 6, c, d} the equation 2 a + 2 fe = 2 C + 2^ /toWs if and only if {a, 6} = {c, ci}. 
Therefore, \So,n\ = and hence K p = 0. 

The next Theorem shows that if k p = 2 P then the limit eigenvalue distribution is the famous 
Marchenko-Pastur distribution. 

Theorem 12. Let k p = 2 P then for every r and p 6 V(r) the coefficient K p = if the partition is 
crossing. Hence 

m r = \NC(r)\ 

the number of non-crossing partitions and p is the Marchenko-Pastur distribution 
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Figure 7. The polytope (x, y, z) e [0, l] 3 such that ft <x + y — z <l. 



Proof. We already observed that K p = 1 iff the partition is non-crossing. Therefore, we need to 
show that for every crossing partition K p = 0. Let p 6 V(r) be a crossing partition with blocks 
{B%, B2, ■ ■ ■ , B\ p \}. Let / be the set of indices such that for i £ I the block Bi does not cross any 
other block Bj. Then we can decompose p as p = p\ U P2 where P2 = Ui^jBi is the union of all 
the non-crossing blocks. Then by the definition of K p we see that K p = K pi K P2 . Now we need to 
show that K pi = 0. Up to a re-enumeration, if necessary, we see that p\ £ V(s) where s < r. By 
definition every block of p\ crosses at least another block. For every n-tuples of positive integers 
(ax,..., a n ) and (bx, ... , b n ) the equation 

2 fl i + 2 a2 + . . . + 2 a " = 2 bl + 2 fe + . . . + 2 bn 

implies that {ax, ■ ■ ■ ,a n } = {bx, ■ ■ ■ , b n }. Hence, every equation in S Pi< n eliminates at least two 
variables and therefore l-Sp^jvl = 0(N s+1 ~ 2 \ pi \). This implies that 

K P1 = lim [%4i , = 
finishing the proof. □ 

In Figure [8j we see the histogram of the matrix VV* for the sequence k p = 2 P and N = 100 
over 1000 trials in comparison with the Marchenko-Pastur distribution. As it can be appreciated 
even for ./V as small as 100 the two are not to far apart. In Figure [9j we see the histogram of 
the eigenvalues of VV* for k p = p — 1 and N = 100. The case k p = p 2 is an interesting one 
(as well as the cases k p = p a ). At the moment we don't understand what is the limit eigenvalue 
distribution for this sequence. For instance, is it true that K p = for every crossing partition? Is 
it true that K p = for the partition p = {{1, 3}, {2, 4}}? In a private communication with Prof. 
Carl Pomerance it was indicated that | ff p jy | is of the order 0(N 2 log(iV)). However, we are not 



providing a proof of this fact. In Figure 10, we show the values of |5 Pi Ar|/iV 3 as a function of N 
and we compare it with the case k p = 2 P . 
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Eigenvalue Distribution of a Generalized Random Vandermonde Matrix 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 8. The blue graph is the histogram of eigenvalues of the matrix VV* for 
the sequence k p = 2 P and N = 100 over 1000 trials. The red curve is the Marchenko- 
Pastur distribution. 

Eigenvalue Distribution of the Classical Random Vandermonde Matrix 

1 i i i i i i 

0.9 

0.8 




Figure 9. Histogram of the eigenvalues of matrix VV* for the sequence k p = p — 1 
and N = 100 over 1000 trials. 
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